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Abstract. We have performed an optical observation 
campaign on PKS 2155-304, whose aim was to determine 
the variability properties of this object on very short time 
scales in several photometric bands. We detected variabil¬ 
ity on time scales as short as 15 min. 

The Fourier properties of the light curves have been 
investigated using structure function analysis. The power 
spectra are well described by a power-law with an index 
—2.4. It is compatible with the index found in the X-ray 
domain. The value of this index shows that the light curves 
cannot be generated by a sum of exponential pulses. Us¬ 
ing historical data, we find that the longest time scale of 
variability in the optical domain lies between 10 and 40 
days. 

We find a strong correlation between flux and spectral 
index, which we interpret as the signature of an underlying 
constant component. As a result we do not find evidence 
of spectral variation for the active nucleus in the optical 
domain. 

A lag has been found between the light curves in dif¬ 
ferent optical bands. The short-wavelength light curves 
lead the long-wavelength ones. The amplitude of the lag 
is about 40 min for a factor 2 in wavelength. 

Our results are compared with predictions from differ¬ 
ent models. None of them can explain naturally the set 
of results obtained with this campaign, but we bring out 
some clues for the origin of the variability. 

Key words: BL Lacertae objects: individual: PKS 2155- 
304 - Galaxies: photometry 


1. Introduction 

Very soon after their discovery BL Lac objects were noted 
for their fast variability. As early as in 1970, Racine ob¬ 
served BL Lac with a time resolution of 15 s and recorded 
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variations of 0.1 mag over a few hours. Since then, many 
BL Lac objects and active nuclei of other types are known 
to vary on time scales shorter than one day in the optical 
domain (e.g., Wagner & Witzel 1995). These objects are 
usually classified as BL Lac objects and blazars, and have 
flat radio spectra. The very fast variability and the radio 
spectrum are often interpreted as the result of synchrotron 
emission from a Doppler-boosted jet. 

PKS 2155-304 is the brightest BL Lac object in the 
optical domain, and is consequently very frequently ob¬ 
served. Variability on time scales of the order of one day 
has been observed by Miller & Carini (1991). In a 4-night 
campaign, they obtained a very smooth light curves, which 
is well eye-fitted by a unique sine function with a 3-day 
period. A drop of 10% has been detected in the ultravi¬ 
olet in a few hours (Edelson et al. 1991). Tagliaferri et al. 
(1991) have studied EXOSAT observations of PKS 2155- 
304, and found that it was dominated by low-frequency 
variability, but detected variability on time scales of the 
order of 500 s. An important campaign has been per¬ 
formed in November 1991 with a very broad spectral 
coverage (Urry et al. 1993; Brinkmann et al. 1994; Cour- 
voisier et al. 1995; Edelson et al. 1995, hereafter 1991-1, 
1991-11, 1991-III, 1991-IV, respectively ). The ultraviolet 
light curve showed some evidence of characteristic time 
scale of ~ 0.7 days (1991-1). During this period, the op¬ 
tical and ultraviolet emissions were well correlated, with¬ 
out evidence of lag (1991-IV). The optical observations 
have been obtained from different ground-based telescopes 
(1991-III), which limited the accuracy of the observations. 
Moreover the sampling was not very dense. As a conse¬ 
quence, the optical variability pattern of PKS 2155-304 
has not been fully explored. 

In this paper we describe an optical campaign on PKS 
2155-304 aiming at determining the properties of the vari¬ 
ability of PKS 2155-304, with an emphasis on very short 
time scales. A limited spectral analysis has also been per¬ 
formed in hve-colour photometry to study whether the 
hardness of the spectrum is related to the brightness of the 
source. We also investigate whether a lag can be detected 
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between the different light curves. On the basis of our 
results, we make some inferences on the emission mecha¬ 
nism of the optical radiation in this object (and possibly 
in other BL Lac objects). A preliminary analysis of the 
data discussed here have been presented in Paltani et al. 
(1996). 

2. Observations 

2.1. The campaign 

All observations were made with the 70-cm Geneva tele¬ 
scope at the European Southern Observatory in La Silla, 
Chile. The telescope is equipped with a CCD camera 
using a thick front-illuminated UV-coated GEC P8603 
416x578 device (Blecha et al. 1990). Besides the usual im¬ 
ager mode, the CCD camera can be operated as a multi¬ 
channel photometer. The Hlters that have been used in 
this campaign are the filters U, B and V from the Geneva 
photometric system (hereafter Ug, Bq and Vq to avoid 
confusion with Johnson’s filters; they are however compa¬ 
rable both in central wavelength and in width), a Hlter 
close to Gunn’s R filter (hereafter Rq) and a filter close 
to Cousins’ I Hlter (hereafter Ic). 

The campaign started on July 26, 1995, and was sup¬ 
posed to last for 3 weeks. However, a bad-weather pe¬ 
riod prevented us to observe after Aug 9. This resulted 
in 15 consecutive nights of observation. We observed PKS 
2155-304 as often as possible, with a rate strongly depen¬ 
dent on the filter. Filters Vq and Rg have been most 
frequently used. As the Ug observations were difficult to 
obtain, only one or two of them were performed each night 
to increase the spectral coverage. The exposure times have 
been estimated a priori to obtain a 1% accuracy on the 
flux (without taking into account the absolute flux cali¬ 
bration), which led us to make 2-min exposures in the Vg 
filter. The pointing and tracking limitations of the tele¬ 
scope reduced the operation period to about 5 hours per 
night, with a gap of about 45 min when PKS 2155-304 
was too close to the zenith. Table 1 summarizes the ob¬ 
servations performed during this campaign. 

2.2. Data reduction 

Only the relevant - user specified - areas of the CCD 
are read out, which enabled us to reduce significantly the 
delay between two observations. Images of repeated ex¬ 
posures are stored in a single structure and reprocessed 
off-line as follows. The raw flux of the object is obtained 
by integration within an aperture of ~ 20 pixels (8 arc- 
seconds). The synthetic aperture is centered using a 2D 
profile fitting. The sky background is determined individ¬ 
ually for each measurement using an average of the flux 
outside the integration zone. The (very rare) cosmic rays 
are detected and removed using the comparison with a 
fitted profile. The fiat-field correction is made on the inte¬ 
grated flux rather than on the raw images (this is equiv¬ 


alent to the filtering of the flat fields). The data acquisi¬ 
tion and image processing are carried out within the stan¬ 
dard Geneva software called INTER (Weber 1993). We 
perform differential photometry with a comparison star, 
which frees us from the problem of atmospheric extinction 
(the object and the comparison star are processed in ex¬ 
actly the same way). The parameters of the comparison 
star, which has been well measured in Geneva photome¬ 
try, are: Ug = 12.036, Bg-Vg = 0.104, Ug - Ug = 1.263, 
Rg — Vg = —0.367 and Ic — Ug = —0.752. Its coordinates 
are 02000 = 21h 59m 02.35s, J 2000 = -30° 10'46.5”. In 
order to take advantage of the differential photometry we 
centred the field so that both the object and the compari¬ 
son star lie in the same frame. Unfortunately, many obser¬ 
vations had to be rejected, because one of the two sources 
had partially disappeared from the field of view during the 
observation, the angular distance between PKS 2155-304 
and the comparison star being close to the size of the field 
of view. 

2.3. Flux calibration 

The absolute flux calibration of the Geneva photomet¬ 
ric system (Rufener & Nicolet 1988) was not intended to 
be used for active galactic nuclei. Spectrophotometric ob¬ 
servations of PKS 2155-304 show that its spectrum is a 
featureless power-law (e.g. 1991-III). To obtain absolute 
flux calibration for PKS 2155-304, we estimate the most 
probable power-law spectrum, given by a least-squares lin- 
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Fig. 1. Mean five-colour spectrum. The empty circles are the 
fluxes from the standard calibration, and the black squares are 
the fluxes corrected to obtain a power-law spectrum. The slope 
of the line is —0.73 
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Table 1. General information on the sampling. n{At < 10 min) and n{At < 5 min) are the number of consecutive observations 
separated by less than 10 min and 5 min respectively. /Norm is the correction factor applied to the flux calibration (see Sect. 2.3). 
The basic statistics of the observations are also given, a is the standard deviation of the light curves in physical units and a is 
the standard deviation in percentage of the mean flux. esF is the uncertainty on the flux obtained from the structure functions 
(see Sect. 4) 


Filter 

Exposure 
time (s) 

Number of 
observations 

n(At < 10 min) 

n(At < 5 min) 

/*Norm 

Mean Flux 
mJy 

(J 

mJy 

a 

% 

esF 

mJy 

Ug 

300 

26 

0 

0 

1.037 

13.3 

2.1 

15.4 

- 

Bg 

180 

110 

49 

0 

0.958 

15.5 

2.3 

14.7 

0.12 

Vg 

120 

262 

188 

86 

0.968 

18.9 

2.7 

14.0 

0.13 

Rg 

60 

224 

144 

70 

1.052 

21.8 

2.7 

12.2 

0.09 

Ic 

80 

130 

62 

8 

0.989 

24.6 

3.1 

12.6 

0.14 


ear regression, from the mean fluxes obtained with the 
standard calibration. The systematic deviations from the 
power-law spectrum are then corrected by the appropri¬ 
ate factor (at most 5%). Fig. 1 shows the mean spectrum 
with the systematic deviations, and the correction factors 
for each filter are given in Table 1. 

2.4- Light curves 

Table 1 gives the basic statistics of the data discussed here. 
Fig. 2 shows the light curves in the five filters in logarithm 
of the flux with an arbitrary normalization. All the light 
curves look very similar. Fig. 3 shows details of 4 nights in 
the Vg and Rq filters. Complicated structures appear dur¬ 
ing most of the nights. These structures are characterized 
by time scales much smaller than one day. A quantitative 
description will be given in Sect. 4. From visual inspection 
of the light curves, it appears that the flux variations are 
roughly simultaneous, as observed in 1991-III and 1991- 
IV. This point will be further analyzed in Sect. 5. 

3. Spectral index variation 

The relative standard deviation of the flux, a, decreases 
with increasing wavelength (cf. Table 1); the probability 
that an uncorrelated population produced a Spearman’s 
correlation coefficient larger or equal to the measured one 
(r = —0.9) is about 4%. This is a first evidence of hard¬ 
ening when PKS 2155-304 becomes brighter. 

We have estimated the spectral index a for each obser¬ 
vation with the Ug filter. The spectral index as a function 
of the Vg flux is shown on Fig. 4. There is a clear evidence 
of hardening when the source gets brighter. A Spearman’s 
test shows that the probability that no correlation exists 
between the flux in the Vg band and the spectral index 
is of the order of 10“^ (r = 0.79). Note that this kind of 
correlation is affected by a bias (Massaro & Trevese 1996). 
Our result is free from this bias, because we have followed 
the recommendations given in the above paper: the Vg 
band is close to the centre of the UG-to-Ic spectral region, 
and the Vg flux has not been used in the determination 
of the spectral index. 


Although the correlation is clear, it disappears com¬ 
pletely (Spearman’s r = —0.04) if one takes only the ob¬ 
servations with a Vg flux larger than 18 mJy. The cause 
of this change of behaviour is discussed in Sect. 6.2. 

On Fig. 4 we have also plotted the data from the 1991 
campaign for comparison (1991-III). It is obvious that 
there is a systematic difference between the two series of 
spectral indices, the spectral indices from the 1991 cam¬ 
paign being larger by a constant about 0.08. This differ¬ 
ence is probably due to the differences in the applied ab¬ 
solute flux calibrations, and in particular to the correction 
that we have chosen to apply. Although two observations 
from the 1991 campaign suggest the existence of the same 
correlations for low fluxes, one observation (around 15.4 
mJy) is clearly not compatible with the present data, and 
another (around 15.1 mJy) is marginally compatible. We 
note that at least two other observations depart from the 
relationship obtained in the present campaign. This is dis¬ 
cussed in Sect. 6. 

4. Amplitude and time scales of the variability 

As we have obtained a very dense sampling, we would like 
to explore the light curves in the Fourier space. However, 
we have to avoid the use of the Fourier transform, since 
it would be strongly contaminated by the window func¬ 
tion, a function fully determined by the sampling. Thus we 
prefer another approach: the first-order structure function 
(hereafter we simply refer to the “structure function”, or 
“SF”). The structure function is a tool commonly used in 
time-series analysis (e.g. Rutman 1978, in clock stability 
analysis). It has been introduced in the field of astronomy 
by Simonetti et al. (1985). The structure function is usu¬ 
ally defined in the case of evenly-sampled discrete time 
series. The structure function Sx(k) of an evenly-sampled 
discrete time series {ti = i ■ At; Xi), i = —oo,..., -foo is 
given by: 

Sx{k) =< [xi+k - Xif >, (1) 

where < fi> is the mean of fi over all i. To approximate 
this relationship in the case where the time series is given 
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Fig. 2. Light curves in the 5 filters in arbitrary units. The light 
curves are plotted in logarithm of the flux with an arbitrary 
normalization. The ticks on the y axis indicate an increase in 
flux by a factor 1.5. The x axis is not to scale 



Fig. 4. Spectral index of PKS 2155-304 as a function of the 
flux in the Vg filter. The empty circles are the observations 
from the 1991 campaign (1991-III) 

hy {tj,Xj), j = 1,... ,n with arbitrary tj, we estimate the 
SF in a bin of width 5 for a lag r using the relation: 

^ ( 2 ) 

’ {i,j) I T — 5l2<tj—ti<T+5l2 
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Fig. 3. Examples of nights where very-short-term variability 
is clearly detected. The filter and the night number are given 



Hours after the first observation 


where N(t,5) is the number of couples [{ti,Xi);{tj,Xj)] 
that satisfy the relationship t — 6/2 < tj — ti < t + S/2. 

In the case of the Fourier transform, the effects of a 
bad sampling are important, but they can be formulated 
in an analytical way, and, in the most fortunate cases, 
properly corrected using a deconvolution. For the SF, we 
do not have such a knowledge about the perturbations 
introduced by the sampling. But because of the particular 
sampling we have obtained, we expect that they will be 
much smaller than in the case of the Fourier transform. 
Indeed, for some values of r, e.g. 0.5, 1.5, 2.5,... days, the 
SF will be completely unknown, because no observation 
can be performed 12 hours after an observation with a 
single ground-based optical telescope. On the other hand, 
for some values of r, e.g. 0.1, 1.0, 2.0,... days, the number 
of couples per bin will be very large, and should lead to 
a very good determination of the SF. Thus, by preferring 
the SF approach, we concentrate the information in some 
parts of the SF, instead of spreading it over the complete 
Fourier domain. 

The SFs of the Bq, Vg, Rg, and Ic light curves are 
shown in Fig. 5. The SF of the Ug light curve has not been 
calculated, because of the very small number of observa¬ 
tions in this filter. Two main features can be observed in 
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Fig. 5. a and b. Structure functions of the Bg, Vg, Rg, and Ig light curves: a SFs for very short r, with 5 = 5 min. b SFs 
for long r, with 5 = 10 min. The solid lines are the mean SFs of the simulated light curves, and the dashed lines show the ±1 a 
limits. The SF of the Ug light curve has not been calculated because of the very poor sampling 


the SFs. The first feature is the presence of a horizontal 
branch at short r. This is due to the fact that the vari¬ 
ability of the time series on very short time scales is dom¬ 
inated by the white noise introduced by the measurement 
errors on the fluxes. As the amplitude of a white noise is 
independent of the lag between the two observations, the 
SF of a pure white noise process is constant, with a value 
equal to twice the variance of the white noise. Thus the 
SF at very short t gives an experimental estimation of the 
uncertainties on the flux, esf- These are given in Table 1. 
The Rg light curve, which has the largest accuracy, shows 
variability on time scales as short as 15 min (0.01 days). 


The second feature is a roughly linear increase of the 
logarithms of the SFs with logr with a slope about 1.3. 
An important property of the SF is that the SF of a time 
series whose Fourier power spectrum follows a power-law 
also follows a power-law. As a first approximation, we as¬ 
sume that the light curves have power-law shaped Fourier 
power spectra. As we cannot correct our SFs from the ef¬ 
fect of the sampling, we use the opposite approach: we 
simulate 100 continuous time series whose Fourier power 
spectra follow a power-law with different indices. To con¬ 
struct the time series, we simply add together sine func¬ 
tions with random phases, with the constraint that the 
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power in each frequency bin decreases as specified. This 
is a discrete approximation of a continuous process, but, 
owing to the very large number of sine functions added 
together, we expect that no serious discrepancy will result 
from this method. We then project the simulated light 
curves onto the sampling obtained for the different filters, 
and calculate their SFs. As a result, we obtain the mean 
SF for our sampling of a light curve that is character¬ 
ized by a power-law-shaped Fourier power spectrum with 
a specified index. 

The power-spectrum index that matches most closely 
the observed SFs of the 4 light curves is —2.4, the indices 
between —2.2 and —2.7 being acceptable (we normalized 
the simulated light curves so that their standard devia¬ 
tion are identical to the one of the observed light curves; 
therefore the SFs of the simulated light curves match the 
one of the real light curve at large r for all indices). We 
can also remark that the features in the SF of the light 
curves (e.g. the relative drop around r ~ 2 days) are not 
at all significant when one considers the dispersion in the 
SFs of the simulated light curves. 

5. Delay between the light curves 

The existence of a delay between the light curves at dif¬ 
ferent wavelengths is an important issue, since it has im¬ 
plications on the emission mechanism. We test this point 
using two methods. 

5.1. Cross-correlation 

The cross-correlation is the most standard way to obtain a 
lag between two time series. We use the interpolated cor¬ 
relation function (ICF) introduced in Gaskell & Peterson 
(1987). This method is more appropriate to our sampling 
than the discrete correlation function (DCF, Edelson & 
Krolik 1988), provided that we calculate the correlations 
for lags smaller than 0.2 days. The reason is that the inter¬ 
polation between two observations made during the same 
night is a very good approximation of the actual flux, be¬ 
cause of the small amplitude of variability for very short 
lags (see Sect. 4). Note that Litchfield et al. (1995) have 
found that the ICF method was more efficient than the 
DCF one when dealing with simulated flaring light curves. 

The cross-correlations of the different light curves with 
the Vg and Rq light curves obtained with the ICF method 
all show a very broad peak that we fitted with a Causs 
profile to find its centre (Table 2 and Fig. 6). The peak 
values are in all cases very close to 1, which indicates an 
excellent correlation. It appears that the shorter wave¬ 
lengths are leading the longer ones. To check whether we 
have really detected a lag, we use the set of 100 simulated 
light curves with an index —2.4 described in Sect. 4, that 
we projected onto the samplings of the real light curves. 
About 10% of the correlations produced a central peak 
that could not be htted by a Causs profile. This produced 



-500 0 500 
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Fig. 6. ICF correlation between the Rg and the Bg light 
curves (solid line). The dashed line is the fit by a Gauss profile. 
Its centre is indicated by a dotted line at t = —44.8 min 


aberrant values for the location of the correlation peak. 
Instead of analyzing each of these correlations to obtain 
a correct value for the lag, we simply discarded all the 
values of the lag outside the range [—0.1; 0.1] days. The 
distribution of the remaining lags is compatible with a 
Gauss distribution. The statistics of the location of the 
peaks are given in Table 2. As the dispersions of the dis¬ 
tributions are about 8 times smaller than the width of 
the domain where we considered the lags as valid, we can 
conclude that the rejection of the values out of the do¬ 
main is justified. The mean lags of the simulations are 
different from 0, but significantly smaller than the dis¬ 
persion, and are thus neglected. Fig. 7 shows the results 
of the cross-correlation analysis. The significances of the 
individual lags are small (< 2a). However two points in¬ 
dicate that the lags are most probably real. First, the lag 
increases monotonously from the Bq to the Ic filters. The 
probability that this happens by chance is the probability 
to obtain a 100% correlation with a Spearman’s test on 4 
points, which is about 2%. Moreover, this has been real¬ 
ized twice independently, since the same observation can 
be made for the lags with the Rq light curve. The second 
point is that the lags between the Bq light curve and the 
Vg and Rg light curves are compatible. This is also true 
for the lags between the Ic light curve and the Vg and 
Rg light curves. 
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Table 2. Results of the cross-correlations of the observed and 
simulated light curves. A negative lag means that the second 
light curve is leading. The last two columns refer to the simu¬ 
lated light curves 


Filters 

Observed lag 
(min) 

Mean lag 
(min) 

Dispersion 

(min) 

Vg-Bg 

-28.2 

+ 14.7 

37.9 

Vg-Rg 

43.5 

- 9.5 

36.9 

Vg-Ic 

58.2 

- 1.4 

32.3 

Rg-Bg 

-44.8 

+ 12.7 

40.4 

Rg“Ic 

23.2 

+ 7.0 

38.2 


Table 3. Results of the global minimization. A negative 
lag means that the second light curve is leading. The “signifi¬ 
cance” is the departure from a zero lag in units of the standard 
deviation for the assumptions = 1 and = 2.5 


Filters Lag Best x^/ Ax^ = Significance 



4000 6000 8000 

Wavelength [ A J 


Fig. 7. Lag between the light curves. The black squares are 
the lags with the Vg light curves (the dotted line is the zero 
lag line). The error bars are the dispersions of the lag in the 
simulations. The black triangles are the analogous of the black 
squares, but with respect to the Rg light curves, except that 
a lag of 43.5 min has been added (the lag between the Vg and 
the Rg light curves). The dashed line is the zero lag line for 
the Rg light curve. The empty circles are the results of the 
global x^ minimization method, and give the lag with the Vg 
light curve. The thick error bars correspond to Ax^ = 1, and 
the thin ones to Ax^ = 2.5 


5.2. Global minimization 

The global minimization method has been introduced 

by Press et al. (1992) to determine the delay between 
two images of a gravitationally-lensed quasar. The princi¬ 
ple of the method is as follows: Assuming two unevenly- 
sampled time series (tf, Xj), z = 1,..., n and (tj, yj), j = 
1,... ,m, we concatenate them to form a new time series 
ZTo,i), i = 1,... ,n + m, one of the two time series 
being shifted in time by an arbitrary lag tq . If the covari¬ 
ance properties of the initial time series are known, one 
can check whether the concatenated time series has the 
same covariance properties. This can be made by calcu¬ 
lating: 

X" = - aEf (C-i) (zr, - aE), (3) 

where C is the total covariance matrix of the process (in¬ 
cluding measurement noise), in analogy with the usual x^ 
definition: x^ = X]r=i The term z — aE mean 

that we subtract a constant a to all components of the 
Zt-o vector (E is the vector unity), a is the value that min¬ 
imizes the X? in the above equation for an assumed lag, 
and is given by: 

ET (C-i) E ^ ’ 

It has the purpose of eliminating the bias introduced by 
the fact the both the mean and the variance measured 
in our data are usually different from their real values. 
This correction even allows us to use the method for low- 
frequency-divergent time series. 

The covariance matrix is constructed on the basis of 
the most complete structure function drawn as a solid line 
in Fig. 5b; it is important to notice that this covariance 
matrix is very well constrained by our data. The details of 
the construction of this matrix are to be found in Rybicki 
& Press (1992) and Press et al. (1992). A drawback of the 
method is that it works only if the Xi and yi data are 
two realizations of the same time series, apart from the 
existence of a lag, and a possible difference in mean value. 
We therefore have to transform the flux in one filter into 
the flux that we would have observed in the other filter. 
We have observed that the flux in a given filter can be 
very well approximated by a constant term plus a linear 
function of the flux in another filter, all the correlation 
coefficients being larger than 0.995. 

The delay is then obtained by minimizing x^('^o)- The 
results are given in Table 3, and are plotted on Fig. 7. We 
see that the lags are essentially consistent with the ones 
found with the ICF method. As it does follow a real x^ 
statistics, we can in principle estimate the uncertainties on 
To by the condition = 1. However Press et al. (1992) 
used the Ay^ = 2.5 after running a set of Monte-Carlo 
simulations. As the computation time is huge, we did not 
try to find the uncertainties on the lags with our simulated 
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light curves. But, contrarily to the situation encountered 
by Press et al. (1992), the values of the lag that we have 
found are in the range where the method should work well; 
therefore we believe that the = 2.5 assumption is 
very conservative. The uncertainties are in any case much 
smaller than in the previous case. The signihcances are 
quite high for all the lags, even if one uses the conservative 
assumption. 

The probability to obtain such values are respec¬ 
tively 12%, 0.8%, and 0.07%. The last two values are very 
improbable. This may be due to the fact that this method 
compares only identical time series. We used experimen¬ 
tal relationships to transform the fluxes in the hlters into 
“equivalent Vq huxes”. This can explain easily the high 
y^ values. 


ist, must take place on time scales even shorter than those 
investigated here, i.e. about 15 min. At least down to this 
limit, variability is produced by a wavelength-independent 
mechanism, for instance geometrical. This possibility will 
be rediscussed below. 

6.1.1. Power-spectrum at very high frequencies 

The structure functions of Fig. 5 are all dominated by the 
white noise introduced by the measurement uncertainties 
for T < 10 min. It means that we can only put an upper 
limit to the minimum variability time scale in this object. 
In addition, the structure function analysis tells us that, if 
the time between two consecutive observations is tq, the 
variability in the Vg band from the hrst observation to 
the second one is: 


6. Discussion 

Emission of BL Lac objects is generally supposed to orig¬ 
inate from synchrotron radiation (at least at low fre¬ 
quency). Other mechanisms for variability have been pro¬ 
posed, like gravitational micro-lensing or geometrical vari¬ 
ations. We are not going to compare detailed models with 
the results obtained here, because of their complexity and 
of the number of free parameters involved. However our re¬ 
sults are strong enough to constrain the qualitative prop¬ 
erties of the models. 


/1.4- 10-4 

V 2 



mJy, 


( 5 ) 


where tq is expressed in min. The accuracy should be of 
the order 0.02 mJy, i.e. 0.1 % of the flux, to detect vari¬ 
ability on time scales of 5 min. The above formula can help 
in the planning of future missions aiming at probing the 
smallest variability time scale in PKS 2155-304. Although 
very-short-term variability does exist in PKS 2155-304, its 
amplitude is very low. This is due to the steepness of the 
power spectrum. 


6.1. Time-series properties of the light curves 

We have obtained a very good description of the Fourier 
power spectra of the optical light curves of PKS 2155- 
304. We can check the compatibility of our result with the 
observations from the 1991 campaign. Fig. 8 shows the 
two SFs of the V (considered identical to Vq) light curve 
from the 1991 campaign and of our Vg light curves. We 
see that both SFs are compatible, apart from the effect 
of the measurement white noise, which has an amplitude 
10 times lower in our data. This comparison may indicate 
that the light curves of PKS 2155-304 are stationary time 
series, or at least that the spectral properties observed in 
the present campaign are not completely peculiar. 

Another comparison can be made with the result of 
Tagliaferri et al. (1991). They found that the Fourier 
power spectrum of EXOSAT observations of PKS 2155- 
304 follows a power-law with an index —2.5 ± 0.2, fully 
compatible with our result. Even after the removal of the 
linear trend, the index (—1.9±0.4) is still compatible with 
our result. This shows that the X-ray and optical emission 
are strongly related, and that they have probably the same 
origin. However, even in this case, many models predict 
that the optical power spectrum will decrease more rapidly 
than the X-ray one, the short-time scale variability be¬ 
ing attenuated at large wavelength, mostly because of the 
longer cooling time (e.g., in synchrotron radiation). This 
is not observed in our campaign. These effects, if they ex- 


6.1.2. Power-spectrum at very low frequencies 

As our campaign lasted only 15 days, it is impossible to 
investigate the existence of time scales longer than this 
duration. However there exists many other observations 
that can be used to determine the complete power spec¬ 
trum. We have used the data from Griffiths et al. (1979). 
The observations are photographic plates that have been 
converted into B magnitudes. Their data are presented in 
two separate sets: a hrst set contains the data from the 
Harvard photometric collection; it spans about 2 500 days 
and contains 66 observations. The second one contains all 
photographic observations, but the data have been aver¬ 
aged by year; the set spans about 90 years and contains 
44 annual averages. The B magnitudes have been trans¬ 
formed into huxes using a 0-mag hux of 4.27 • 10® mJy. To 
transform the variations into variations in the V band, we 
transformed the B huxes into Vg huxes using the linear 
relationship that we have found between the Bg and Vg 
huxes. 

The SFs of the Harvard data set and of the complete 
photographic data set are also plotted. The uncertainties 
on theses SFs are large, because of the poorness of the 
sampling, but we see that the SFs of the two long-term 
data sets do not follow the increase seen in the present 
campaign. As the complete data set is averaged by year, 
one must consider its SF as a underestimation of the real 
SF, because such an average has great chances of reducing 
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Fig. 8. Structure functions of the light curves in the V 
band from our campaign (black squares), the 1991 campaign 
(crosses), the Harvard data set (empty squares), and the com¬ 
plete photographic data set (black circles). The solid line is a 
power-law with an index 1.3 with a break around 20 days and 
a contamination with a white noise with a standard deviation 
of 0.13 mjy. The dashed line is the same power-law, but con¬ 
taminated by a white noise with a standard deviation of 1.34 
mJy. This clearly shows the increase of the accuracy reached 
in our campaign 

the variability. Taking this point into account, this SF is 
compatible with the SF of the Harvard data. The tentative 
SF drawn on Fig. 8 indicates that the longest variability 
time scales in PKS 2155-304 lies between 10 and 40 days, 
apart perhaps from time scales longer than 100 years. 

6.1.3. Origin of the power spectrum 

The variability of PKS 2155-304 is small, and the light 
curves are very different from what is observed in, for in¬ 
stance, OJ 287 (e.g. Sillanpaa et al. 1996), which has an 
amplitude of variability of at least 4 mag, i.e. a factor 40. 
Actually the variability of OJ 287 appears in sudden bright 
“flares”, which could show that the mechanisms that pro¬ 
duce the essential of the variability are very different from 
those at work in PKS 2155-304 (e.g. the binary black hole 
model Sillanpaa et al. 1996). Apart from the distinction 
between “flaring” and “non-flaring” objects, we do not 
know whether the Fourier properties of PKS 2155-304 
are typical of BL Lac objects. We have indeed here the 
most extensive (in term of spectral coverage) study of the 
power spectrum of any active galactic nucleus. The most 
important point is that all the data are compatible with 


a power-law-shaped power spectrum with an index —2.4, 
and and low-frequency cut-off at frequencies as high as 
(10-40 days)-b 

Many other time series in the physical world have 
a power-law-shaped Fourier spectrum (e.g. in electronic 
systems, or, most surprisingly, in the first Brandenburg 
Concerto from J.S. Bach, Voss & Clarke 1975), but no 
convincing physical mechanism has ever been proposed 
(Press 1978). Several physical processes have been pro¬ 
posed for X-ray light curves of Seyfert galaxies. The un¬ 
derlying idea of most of them can be formulated by the 
“mechanical model” of Halford (1968). He has shown an¬ 
alytically that a superimposition of any “reasonable” (as 
defined in his paper) time-dependent perturbations with 
different time scales and amplitude can generate any spec¬ 
tral index in the Fourier domain. As our index is smaller 
than —2, it means that we can constrain the properties 
of the perturbations. It requires indeed that the Fourier 
power spectrum of the perturbations decreases with an 
index at most —2.4. Therefore exponential pulses, whose 
Fourier power spectrum decreases with an index —2, are 
excluded. A candidate for the events could be the injec¬ 
tion of packets of relativistic electrons, which cool by emit¬ 
ting synchrotron radiation. Variations in the parameters 
of the packets (size, energy distribution of the electrons, 
...) could produce the diversity of amplitudes and time 
scales required to obtain a power-law-shaped Fourier spec¬ 
trum (but see Sect. 6.2 for the problem of spectral index 
variation). In this case the low-frequency cut-off is related 
to the maximum cooling time of the packets, if the “birth 
times” of the packets are completely independent from 
each other. 

Another possibility of producing a power-law-shaped 
Fourier spectrum is by summing many periodic functions 
with different periods. Following the idea of Camenzind 
& Krockenberger (1992), a “knot” can have a helical mo¬ 
tion around the magnetic field, which produces periodic 
flares, because of the periodic variation of the Doppler am¬ 
plification towards the observer. A large enough number 
of knots with different gyration radii could fill the power 
spectrum to produce a power-law. Camenzind & Krocken¬ 
berger (1992) found that the typical time scales induced 
around a 10® Mq black hole should be of the order of 1 
day. However it seems plausible that smaller time scales 
can be obtained with modifications of the geometry of the 
system, like the angle between the jet and the observer, 
the velocity of the knot, or the mass of the black hole. 

We can further add that it seems rather improbable to 
us that gravitational micro-lensing can explain the light 
curves observed in our campaign. Even though Kayser 
et al. (1989) have shown that the light curves generated 
by a micro-lensing foreground galaxy can be very com¬ 
plex and have broad Fourier power spectra, their shapes 
show successions of sharp flares, which do not appear in 
the light curves from this campaign or from the 1991 cam¬ 
paign (1991-1, 1991-H, 1991-HI). 
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Because variability is essentially geometrical (provided 
that the cooling times of the packets are significantly 
longer that the time scale of the geometrical variability), 
these last two models predict a power spectrum mostly 
independent of the frequency, as observed. 


6.2. Spectral behaviour of PKS 2155-304 

A correlation between spectral index and flux is clearly 
seen in our data. It has however not always been the case 
in other studies on PKS 2155-304. CCD optical observa¬ 
tions from Zhang & Xie (1996) indicate that there is no 
correlation between the B-V colour index and the V mag¬ 
nitude. The complete archive of lUE spectra shows that 
the variability increases when the wavelength decreases 
(Paltani & Courvoisier 1994), but Edelson (1992) found 
using a large part of this archive that the flux and the 
spectral index were not correlated. On the other hand he 
found a (not very clear) correlation in MRK 421. In the 
1991 campaign the constant spectral index hypothesis was 
favoured in the optical data (1991-III), while a hardening 
of the spectrum when the source brightens was the gen¬ 
eral behaviour observed in the ultraviolet domain (Urry 
et al. 1993), but the details were very complex and do not 
support the clear correlation found in our data. In other 
BL Lac objects, the correlation has been found to be very 
significant (e.g., OJ 287: Gear et al. 1986; AO 0235-1-164, 
PKS 0735-bl78, 1308-b326: Brown et al. 1989). On the 
other hand, in a study of 6 BL Lac objects, Massaro et al. 
(1995) found only two cases of positive correlations. 

Examining Fig. 4, it appears that the spectral index is 
compatible with a constant as soon as the flux in the Vq 
band is larger than 18 mJy. An explanation for the drop 
of spectral index at low fluxes could be that the fluxes are 
contaminated by a constant emission, e.g. from the host 
galaxy, that is fainter and steeper than the BL Lac com¬ 
ponent. PKS 2155-304 was rather weak during this cam¬ 
paign. Indeed the mean V flux obtained here is about 15 % 
smaller than the one obtained during the 1991 campaign. 
Moreover the V flux was below 18 mJy during about the 
half of the campaign, while it was below this value dur¬ 
ing about only 20 % of the 1991 campaign. This, together 
with the much better signal-to-noise ratio obtained in our 
campaign, could explain why the drop of the spectral in¬ 
dex at small flux has not been observed during the 1991 
campaign. A further argument in favour of the existence 
of an underlying component comes from the remark made 
in Sect. 5.2: The best way that we have found to transform 
a flux into a Vq flux is through a linear relationship with 
a constant that increases with the wavelength] a compara¬ 
ble observation in Seyfert galaxies has been interpreted by 
Paltani & Walter (1996) as the signature of an underlying 
component. Therefore we cannot exclude the possibility 
of a constant (or weakly varying) spectral index, and it is 
the interpretation that we favour. 


A few optical observations from the 1991 campaign 
are in disagreement with the above interpretation that the 
variations are achromatic (apart from the existence of an 
underlying component). It is also the case for several of the 
campaigns on PKS 2155-304 cited above. However in the 
two optical campaigns, most of the observations are com¬ 
patible with the constant spectral index hypothesis (note 
that a drop of the colour index can be marginally observed 
in Zhang & Xie (1996)). It suggests that the variability is 
mostly achromatic (or weakly chromatic), with some ca¬ 
sual significant excursions of the spectral index from the 
average value. 

Achromatic variability cannot be easily accounted for 
in terms of synchrotron radiation, because it would require 
a very fine tuning between injection terms and electron 
losses. On the other hand small variations of the angle be¬ 
tween the jet and the observer, or of the bulk Lorenz factor 
can generate achromatic variability. Gravitational micro- 
lensing by the stars of a foreground galaxy has already 
been invoked to explain the variability of BL Lac objects 
(Schneider & Weiss 1987; Kayser et al. 1989). However, as 
mentioned above, variability is not really achromatic, but 
rather the spectral index has little variations, which are 
uncorrelated with the flux. This cannot be easily explained 
either by geometrical variations alone or by micro-lensing. 
As the flux variability is rather small, it may be that spec¬ 
tral index variation is small, if a large number of packets of 
electrons with different energy distributions participate in 
the optical emission. Moreover a campaign performed in 
May 1994 that included ASCA observations showed that 
the variability in the hard X-ray domain (above the range 
covered by ROSAT in 1991-H) was quite different both in 
amplitude and in time scales from what has been observed 
with lUE (Urry 1996). 

6.3. Delay between the light curves 

In 1991-IV, Edelson et al. already observed a lag of 2- 
3 hours between the ROSAT X-ray light curve and the 
lUE 1400 A light curve. We confirm the existence of lags 
between the different wavelengths, in the sense that the 
light curves at shorter wavelengths lead the light curves at 
longer wavelengths. The delays between the optical light 
curves are about 40 min for a factor 2 in frequency. This 
value is actually very close (and compatible, owing to the 
large uncertainties on these values) to the amplitude of 
the lag between the ROSAT X-ray light curve and the 
lUE 1400 A light curve (a factor ~ 50 in frequency). 

Inhomogeneous jets, i.e. where electrons are injected at 
the base of the jet and emit most of the synchrotron ra¬ 
diation at frequencies decreasing with the distance to the 
base, naturally produce a lag in the sense observed here 
(but see Sects 6.2 and 7). It has however been noted in 
1991-IV that the amplitude of the X-ray-ultraviolet lag is 
inconsistent with models where electrons are injected and 
not reaccelerated, because the life time of a X-ray-emitting 
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electron is extremely short (of the order of the second in 
the model of Ghisellini et al. (1985)). The problem also 
exists in the optical domain, where the lag is comparable 
to the one between the X-ray and the ultraviolet emis¬ 
sions, and the life time of the optical-emitting electrons 
has only increased by a factor of the order of 10. 

This problem can be solved if the lag is due to a bub¬ 
ble of electrons that becomes progressively optically thin 
at longer wavelengths (Marscher & Gear 1985). We note 
however that we disagree with the observation made in 
1991-IV that the lag is shorter than the shortest variability 
time scales. The 2-3 hours observed in the X-ray domain 
are at least 8 times longer that the shortest variability 
time scale observed in the Rg light curves, and assuredly 
only the measurement noise has prevented us to detect 
even shorter time scales. This is difficult to explain with 
this model. 

Although stated in the literature (e.g., 1991-IV), we 
do not consider plausible that gravitational micro-lensing 
models can account for the existence of a lag. The mech¬ 
anism proposed in 1991-IV should produce a broad, but 
symmetrical with respect to r = 0, correlation peak. 

Note that the presence of lag between two light curves 
automatically implies spectral-index variation. In PKS 
2155-304, the short-time scale variability is sufficiently 
small, so that the variation of the spectral index due to 
the lag remains undetected. 

7. Summary and conclusion 

7.1. Results on the variability of PKS 2155-304 

We have performed several types of analyses to extract 
the physical properties of the optical light curves of PKS 
2155-304. The main results are the followings: 

1. The structure function of the optical variability is com¬ 
patible with the one expected for a power-law-shaped 
Fourier power spectrum with an index —2.4. 

2. The shortest variability time scale observed is about 
15 min, but it is an upper limit, as the variability is 
dominated by the measurement white noise below this 
time scale. 

3. The optical spectrum hardens when the source bright¬ 
ens. We conclude however that this is probably due to 
an underlying stable component. As a consequence the 
variability of the BL Lac component in PKS 2155-304 
is considered to be achromatic during this campaign. 

4. The fluxes at different wavelengths are very well cor¬ 
related. 

5. There is a lag between the different light curves. The 
light curves at longer wavelengths are delayed with re¬ 
spect to the light curves at shorter wavelengths. The 
amplitude of the delay is about 40 min for a factor 2 
in frequency. 

To these results, we can add those obtained from the 
comparison with other data: 


6. The maximum variability time scale is comprised be¬ 
tween 10 and 40 days. There is no evidence of variabil¬ 
ity on time scales between 40 days and 100 years. 

7. The shapes of the power spectra in the optical and 
X-ray domain are compatible. 

8. The optical spectral index has been observed to vary 
slightly without correlation with the flux during other 
campaigns. 

9. The lag between the X-ray and the ultraviolet light 
curves is 2-3 hours for a factor 50 in frequency, i.e. 20- 
30 min for a factor 2 in frequency, which is very close 
to the lag in the optical domain. 

7.2. General conclusions 

We have considered the global properties of several models 
that could explain the variability of PKS 2155-304: syn¬ 
chrotron radiation from an inhomogeneous jet (Ghisellini 
et al. 1985), synchrotron radiation from an shocked bubble 
(Marscher & Gear 1985), variability induced by geomet¬ 
rical variations (Camenzind & Krockenberger 1992), and 
variability induced by gravitational micro-lensing events 
(Schneider & Weiss 1987; Kayser et al. 1989). The follow¬ 
ing conclusions can be drawn from the above results. They 
are sufficiently general to be independent of the details of 
the models. 

Power spectrum origin Points (1), (2), and (6) give a very 
complete knowledge of the variability properties of the 
optical light curve. The meaning of the low- and high- 
frequency limits of the power spectrum cannot be under¬ 
stood without a knowledge of the variability mechanism. 
This is also true for the shape of the power spectrum. But 
we can exclude models where the light curve is generated 
by the sum of a large number of exponential pulses, while 
other shapes of pulses are fully viable. We can moreover 
note that some models are less susceptible of generating 
the kind of light curves and power spectra observed here, 
in particular this seems to be the case for the micro-lensing 
model. 

Power spectrum constancy Points (4), and (7) indicate 
that, if the optical power spectrum got through a high- 
frequency filter, the low-frequency cut-off of this filter 
would be at very high frequencies (> (0.01 days)“^). 
We prefer another explanation, namely that the high- 
frequency variability has a geometrical origin. 

Gravitational micro-lensing models can be ruled out by 
Points (3), (4), and (5): Point (3) indicates that the size 
of the blue-light emitting region is the same that the one 
of the near-infrared-light emitting region, but Point (5) 
shows that these regions cannot be identical, and Point (4) 
would imply that all lenses cross first the blue-light emit¬ 
ting region, and about one hour later the near-infrared- 
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light emitting region with the same amplihcation param¬ 
eters. This is, of course, extremely improbable. 

Amplitude of the lag Points (5), and (9) show that the 
lag observed between the light curves is much larger than 
the one predicted by inhomogeneous jet models without 
reacceleration. A model in which a dense bubble subject 
to synchrotron self-absorption becomes progressively op¬ 
tically thin can however explain this lag. 

Lag - variability time scale problem We must remark that 
Points (3), (4), and (5) are difficult to reconcile for any 
model. Models where variability is “geometrical” can ac¬ 
count for the achromaticity of the optical variability, but it 
appears difficult to explain a lag signihcantly longer than 
the minimum variability time scale. This favours an inter¬ 
pretation where low- and high-frequency variabilities have 
different origins (e.g. synchrotron cooling and geometrical 
variations respectively). 

(Quasi-) achromaticity of the light curves Synchrotron 
emission does not naturally produce achromatic varia¬ 
tions, as observed in Point (3). If synchrotron radiation is 
the mechanism that produces the optical emission of PKS 
2155-304, the slight spectral index variations (Point (8)) 
make us think that the optical emission is produced by a 
number of bubbles sufficiently large to “dilute” efficiently 
the spectral index variations induced by each bubble. 

The variability picture of PKS 2155-304 remains a very 
difficult puzzle, but the results given in this paper are 
very constraining. To our mind, the most interesting issue 
raised here is to find a mechanism that can explain simul¬ 
taneously important lags and well-correlated variability 
on very short time scales. But we keep also in mind that 
other events not described here, because not observed dur¬ 
ing this campaign, exist in PKS 2155-304, like the ASCA 
flare (Urry 1996). 
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